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Abstract 

The market events of 2007-2009 have reinvigorated the search for re- 
alistic return models that capture greater likelihoods of extreme move- 
ments. In this paper we model the medium-term log-return dynamics in 
a market with both fundamental and technical traders. This is based on 
a Poisson trade arrival model with variable size orders. With simplifica- 
tions we are led to a hybrid SDE mixing both arithmetic and geometric 
Brownian motions, whose solution is given by a class of integrals of expo- 
nentials of one Brownian motion against another, in forms considered by 
Yor and collaborators. The reduction of the hybrid SDE to a single Brow- 
nian motion leads to an SDE of the form considered by Nagahara, which 
is a type of "Pearson diffusion", or equivalently a hyperbolic OU SDE. 
Various dynamics and equilibria are possible depending on the balance 
of trades. Under mean-reverting circumstances we arrive naturally at an 
equilibrium fat-tailed return distribution with a Student or Pearson Type 
IV form. Under less restrictive assumptions richer dynamics are possi- 
ble, including bimodal structures. The phenomenon of variance explosion 
is identified that gives rise to much larger price movements that might 
have a priori been expected, so that "25a" events are significantly more 
probable. We exhibit simple example solutions of the Fokker-Planck equa- 
tion that shows how such variance explosion can hide beneath a standard 
Gaussian facade. These are elementary members of an extended class of 
distributions with a rich and varied structure, capable of describing a wide 
range of market behaviours. Several approaches to the density function 
are possible, and an example of the computation of a hyperbolic VaR is 
given. The model also suggests generalizations of the Bougerol identity. 

Keywords: Market microstructure, fundamental trader, technical trader, Stu- 
dent distribution, t-distribution, Skew-Student, Pearson Type IV, Fokker-Planck 
equation, stochastic differential equation, partial differential equation, credit 
crunch, variance explosion, Bougerol identity, Asian options, exponentials of 
Brownian motion. 
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1 Introduction 

'Technical analysis is an anathema to the academic world. ' 
Burton Malkiel, in "A Random Walk Down Wall Street" . 

The fat-tailed and non-normal distribution of share-price returns has been 
well known for decades. Mandelbrot [26] and Fama [18] noted the excess kurtotic 
nature of equity returns in the 1960s. The need for non- (log) Gaussian has been 
widely recognized by numerous authors, and the events of 2007-2009 have made 
it very clear that extreme share-price movements arc much more likely than in a 
simple log-normal model associated with geometric Brownian motion. Despite 
the fat-tailed behaviour being known for 40 years many market participants 
cling to the Gaussian or near Gaussian picture, and attempt to explain away the 
large movements while remaining within that picture. The CFO of Goldmans, 
in a Financial Times article, famously attempted to excuse the implosion of 
Goldman's hedge funds with the comment 1 

"We were seeing things that were 25-standard deviation moves, several days 
in a row". 

These were not events that can be written off as un-modcllablc due to their 
extreme unlikelihood - people were just thinking of the wrong, Gaussian, distri- 
bution based on an unrealistic world- view, while simultaneously discounting the 
possibility of both major shifts in the mean and sudden increases in variance. 
For example, the likelihood of a single 25cr or worse event, without allowing for 
a major shift in the distributional mean or variance, is 

• about 6 x 10~ 138 in a Gaussian picture; 

• about 4x 10~ 6 in a Student-t picture (four degrees of freedom, as estimated 
for global indices in [19]); 

• about 1/625 for a toss of a very unfair coin where the probability of a 
head is 1/625. 2 

The estimation of the observed return distribution is a matter for proper 
statistical analysis. The emergence of some form of fat-tailed or indeed power- 
law decay in the tails is a robust feature on various time-scales. The leading 
contenders for the best fit are often a variance-gamma (VG) distribution or a 
Student t, with some evidence for normal-inverse Gaussian. A useful approach 
is to do a maximum likelihood estimate within the generalized hyperbolic family. 
See, for example, the recent study by Taylor et al on South African index data 
[41], Fergusson and Platen [19], and references therein. Extreme events are 
much more common in such models. The work of [19] found a Student with 
about four degrees of freedom (y = 4), and the work of [41] ranked a Student as 
second most likely with higher values of v. General non-integer low values of v 
may well be of interest in financial analysis for short time- scales. Work cited in 
[22] suggests that very short term returns exhibit power law decay in the PDF. 
The values of v reported in [22] take values in the range 2 to 6. So this leads us 
to consider not only small integer values of v : 2 < v < 6 but also non-integer 
v. In normal statistical estimation use the v parameter is usually assumed to 
be integer as it relates to sample size, but there is no underlying mathematical 
reason for this restriction. 

The huge losses in the markets in the fallout of the credit crunch were due 
to several factors, which include not only the very real loss in value of financial 

1 David Viniar, Goldman's chief financial officer, as quoted by Peter Larsen, Financial 
Times August f3 2007. 

2 This example is based, with thanks, on an anonymous observation at 
http: //worldbeta.blogspot . com/2007/08/really- with- seth- and- amy-part-ii .html 
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institutions who might, for example have substantial CDS exposure or exposure 
to genuine economic downturn, but also to the fear of traders panicking to 
offload falling assets. 

But long before the recent panic, participants in financial markets traded 
according to diverse strategies. Two common approaches are the so-called "fun- 
damental" and "technical" trading strategies. Traders in the first group are 
largely interested in measure of value and traders in the second group are inter- 
ested more in the dynamics of the price (and volume) history. In this note the 
impact of having these two types of trades is analyzed employing a simple trade 
arrival model with Poisson characteristics, with variable size orders. In general 
we obtain a large class of price evolution models. 

The approach presented here has conceptual links with other approaches. 
In particular, Brody, Hughston and Macrina ('BHM') [10] have introduced the 
notion of information-based asset pricing, paying rather more attention to the 
filtration component of processes than has hitherto been developed. In the 
model presented here we package information received since the start of a trad- 
ing period into the price history since the trading start, on the basis that mar- 
ket participants who believe in technical trading are at work. 3 Our linearized 
model results in an SDE that is, at least at the conceptual level, a linear price- 
information correction to standard Brownian motion of log-returns. It remains 
to be seen whether this notion can be given a rigorous mathematical basis within 
the BHM framework. 

We consider the dynamics on a medium time scale which may be consid- 
ered as assessing the distribution of returns on a daily basis as the aggregate 
of many intra-day trades. Under certain simplifying assumptions, correspond- 
ing to linearization of price impact functions, linearization of technical trading 
criteria, and a "many trade" limit, we are lead to an SDE that is a hybrid of 
arithmetic and geometric Brownian motion. Various dynamics and equilibria 
arise from this SDE depending on the balance of trades. Under mean-reverting 
circumstances we are lead naturally to an equilibrium fat-tailed return distribu- 
tion with a Student t or skew-Student form, with the latter defined within the 
framework of "Pearson diffusions" defined by Forman and S0rensen [20]. The 
construction of the standard Student distribution from an SDE has also been 
recently considered by Steinbrcchcr and Shaw [39] . The model is capable of still 
richer dynamics, and in more generality leads to a hyperbolic extension of the 
OU process. In general form a non-linear diffusion with or without jumps is pos- 
sible. One effect of note is the simultaneous explosion in the variance coupled 
to the emergence of a non-Gaussian distribution, and such distributions hiding 
their character under a Gaussian mask. When there is significant momentum 
bimodal densities may arise. 

This note is not to first to propose a model with Student equilibria. Indeed 
the work by Nagahara [31] , that postulated a certain Pearson diffusion, pre-dates 
our analysis significantly, and the associated density function was given by Wong 
[43] many years ago (it is Wong's "Type 2E"). The work by Carmona, Petit and 
Yor [11], set in the context of Levy processes, contains closely related structures. 
However, the full financial market and risk significance of this model, and its 
potential justification in terms of price feedback via technical or similar trading, 
does not appear to have significant penetration in the financial mathematics 
community, and the technology associated with managing the resulting densities 
is not developed. In this current note we provide: 

3 This of course makes no judgement as to the wisdom, mathematically, statistically or 
otherwise, of technical trading, but does acknowledge its presence. In fact all we really need is 
the presence of price-sensitive orders. In reality the weight of technical trading is growing, with 
complex program trades taking up an increasing fraction of the market. Even fundamental 
traders will demand efficient execution and hence acquire a degree of price-sensitivity. The 
model developed here is a simple, linear and Markov representation of this market complexity. 
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1. a sketch of a micro-structure justification; 

2. the connection with the theory of integrals of exponentials of Brownian 
motion; 

3. an alternative density representation, and some simplified forms of the 
density; 

4. an analysis of market states and pathological cases; 

5. an initial approach to the VaR analysis. 

In this way we hope the reader will be persuaded of the rationale, underlying 
simplicity and practical applicability of this model, which produces fat-tailed 
behavior consistent with market behaviour in a simple manner. There is of 
course a jump-extended version of the model incorporating technical trading. 
Here we do not allow for individual "large" trades. Research into more tractable 
forms of Wong's distribution [43, 31] is also needed. We will offer a simple single- 
term formula for the transform here. The theory of integrals of exponentials of 
Brownian motion, and of Levy processes, has attracted much study over the 
years and is of considerable interest from the point of view of both pure analysis 
and practical applications, initially to the Asian option. The reader is referred 
to the articles by Dufresne [16, 15] and to the classic collection of papers by M. 
Yor [44] and the references therein for the detailed mathematical background 4 . 
Several properties of integrals of GBM are also reviewed in [14]. A more recent 
survey of some key ideas is given by Matsumoto and Yor [28], see also [27]. 
For the PDE enthusiast we refer the reader to Dewynne and Shaw [13], where 
an unusually quick derivation of the Asian Call price is given together with an 
volatility series, cf. Zhang [46]. There are many other papers on approximations, 
bounds, trading, discrete averaging but the key papers above and references 
therein are most relevant to the exact representation in terms of integrals of 
exponentials of Brownian motion. For practitioner applications many discrete 
and continuously averaged cases of Asian options may be computed by the 
approach of Vecer [42] . 

In general the probability density functions arising from the model are com- 
plicated. The direct Laplace transform approach presented here offers an alter- 
native representation to that given by Wong [43], and has the advantage that 
several special cases are closed-form computable for both mean-reverting and 
momentum-dominated states with integer drift. This includes both some cases 
(mean-reverting) identified by Wong, as well as simple momentum-dominated 
cases intimately related to the Bougerol identity in both the original form, [8], 
and the beautiful bi- modal variation found by Alili, Dufresne and Yor [2]. It 
is hoped that the representations given here may lead to further results along 
these lines. 

1.1 Analogies from the physical and biological world 

We will work with an SDE that is a hybrid of arithmetical and geometric Brow- 
nian motions. This SDE has many interesting features, and has an interesting 
history in the physics literature 5 . A very closely related analysis has been given 
quite recently in the plasma physics literature [38]. However, the history of 
hybrids of both arithmetical and geometric Brownian motions dates back to the 
1979 paper by Schenzlc and Brand [34]. Also, discrete models in the form of a 

4 The author is indebted to Prof M Yor for pointing out many of the relevant papers on 
related topics - any omissions that remain are due to me. 
5 I am grateful for Dr. G. Steinbrccher for these insights. 
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linear Langevin equation: 



x(t + l) = b(t)x(t) + f(t) 



(1) 



where b and / are both random, were considered in a thread of research originat- 
ing in 1997 with Takayasu et al, who established conditions for the realization 
of power law behaviour [40]. 

On the biological side, May et al [29] have observed the "common ground in 
analysing financial systems and ecosystems" , and identify feedback mechanisms 
as one mechanism for causing catastrophic changes in the overall state of a 
system, with hidden linkages being another cause. Here we are focusing on the 
mechanism of price feedback within the financial context, while adapting the 
familiar financial mathematics framework of stochastic differential equations to 
cope with it. The importance of feedback has also been emphasized by Kambhu 
et al [24], in a major survey encouraging the financial community to draw on 
the lessons from ecology, atmospheric science and other complex systems. It is 
hoped that this current paper will contribute to that discussion. 



The plan of this work is as follows. In Section two the underlying dynamics 
of a market containing both fundamental and technical traders is considered. 
In Section three the model is simplified and linearized, and approximated to 
a tractable model. In general a model with a combination of jumps and pure 
Brownian motion is possible, but we focus on the non-jump case for further 
analysis. In Section four the possibility of a Student t equilibrium is identified. 
In Section five the notion of a "hyperbolic OU" process is established, and in 
Section six the possibility of Pearson Type IV is considered. The fully dynamical 
case is analyzed in Section seven. In particular the phenomenon of "variance 
explosion" is identified. We also demonstrate the nature of the full distribution 
by reference to an explicitly solvable special case, linked to Bougerol's identity. 
In Section 8 we analyze the general case, including examples of both momentum- 
dominated and mean-reversion-dominated markets. This section also explores 
other closed-form distributional solutions and gives a Laplace transform solution 
for a large class of cases, revealing the link to the Legcndrc equation well known 
in physics, and to the Bougerol identity and its bimodal generalization. Section 
nine offers a preliminary market classification based on this approach. Section 
ten explores the VaR implications and Section eleven offers conclusions and 
speculations. 

2 Two types of market participant 

This section is intended to be a sketch of a series of links from basic trading 
ideas through to a stochastic differential equation, and is intended in that spirit 
only. There are aspects where a more rigorous mathematical treatment will 
ultimately be needed, such as the details of the SDE limiting process. However, 
we think it is important to provide such a sketch. The later analysis could 
proceed basely solely on the postulation of the resulting hybrid SDE, but it is 
intended that this paper provide an end-to-end discussion, starting with ideas 
about trades and ending with a density function. 

Consider a market containing an asset with share price St at time t. We 
consider a trading period t € [0, T]. The log-return, x t , on the asset is given by 



1.2 Plan of article 




(2) 
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and we shall be concerned with time intervals over which St ~ Sq(1 + Xt) so 
that no distinction is made between log and linear returns. 

Our approach is to consider that the agents trading in this asset have two 
different types of motivation. 

One set of agents comprise those acting independently of the current value 
of x t . These will include, but will not necessarily be limited, to those trading 
on fundamentals, such as the dividend or earnings yield based on the price Sq. 
Other such traders may be engaged in a portfolio rebalance, or re-hedging a 
derivative position. 

The second set of agents trade on the basis of the value of x t . These are 
our technical traders relative to the time period under consideration. Technical 
trading is based on a number of different algorithms, which may be classified 
according to a "Rumsfeld" scheme: price momentum, mean reversion (known 
knowns) through unknown unknowns such as black box unpublished models 
being run by covert trading operations. We do not know the totality of such 
agents but we do know that they care about x t , or possibly about the history 
{x s \0 < s < t}. In this paper I shall simplify and consider traditional momentum 
and mean-reversion traders who act on the value of x t . There may also be longer 
term technical traders focusing on the ratio of So to some much older price, who 
happen to be at work in the given time interval, but for our purposes they will 
not count as technical as they are not reacting specifically to x t . There may be 
other participants who would not wish to be considered as technical traders, e.g 
those carefully achieving a position by a sequence of trades, but to the extent 
that their willingness to trade depends on x t they are technical. Neither shall we 
enter the debate as to whether technical trading makes any sense, whether price 
increments are independent or whether there are material serial correlations etc. 
etc. It will suffice that people who believe in technical trading are trading. In 
this sense Malkiel's famous objection is irrelevant. 

2.1 Fundamental "buy" orders 

Let us now consider a time interval (t, t + At) C [0, T] with At « T, and that 
orders may be effected in lots of size L. We consider first buy orders based on 
fundamental trading. 

Let Y be the integer-valued random variable denoting the number of such 
trades arriving in time At, and let Ni, i = 1, . . . Y, be the integer-valued random 
variable denoting the number of lots in each buy order. The number of shares, 
Mb, in the total collection of buy orders is M = L x Z, where the random 
variable Z is given by 

Y 

Z = Y,Ni. (3) 

i=0 

We assume that the number of trades is independent of the size of each trade, 
and that the trade sizes are independent and identically distributed 6 . Then the 
probability generating function (PGF) fz(s) of Z is related to the PGFs of Y 
and ./V by 

fz(s) = MM*)) , (4) 
from which elementary PGF theory tells us that 

E\Z\ = E[Y]E[N] = E[Y]n , (5) 

"These are assumptions and ones that might reasonably be questioned. One might consider, 
for example, that in a significant market downturn, there is correlation between having larger 
trades and having more trades. But even in this case we can imagine many small investors 
selling as well. Our goal is to get a model of price feedback, and the elegant compositional 
relationship for PGFs that these assumptions enable allow us to proceed more easily, if not 
in complete generality. 
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where E[N] — n, and 

Var[Z] = Var[Y]n 2 + E[Y]Var[N] . (6) 

We shall now assume that the buy variable Y either follows a Poisson process 
with arrival rate X B , or, more loosely, that the process is sufficiently Poisson-likc 
that we can write 

E[Y] = Var[Y] = X B At . (7) 

In fact what we really need is only that this last equation holds. Y following a 
Poisson process is sufficient but not necessary. So then we have 

E[Z] = X B AtE[N] = X B Atn , (8) 

Var[Z] = \ B At(E[N} 2 + Var[N}) = X B At(n 2 + Var[N}) = \ B AtE[N 2 } . (9) 
It follows that 

E[M B ] = LX B Atn , (10) 
and that the standard deviation of M B , sd(M B ), is 

sd(M B ) = Ly/X B AtE[N 2 } . (11) 

2.2 Fundamental "sell" orders 

This proceeds in the same way. With similar assumptions, including that the 
variable N is similarly distributed, we end up with mean of the number of sales 
in the sell orders as 

E[M S ] = LXsAtn , (12) 
and that the standard deviation of Ms, sd(M$), is 

sd(M s ) = Ly/X s AtE[N 2 } . (13) 

2.3 Aggregation of fundamental trades 

We define 

Mp = M B - M s (14) 
as the net buy volume. The random variable Mp has mean 

E[M F ] = L{X B -X s )Atn . (15) 

and its variance depends on the correlation between arrival rates of buy and 
sell trades. If we assume that fundamental buyers and sellers are acting on 
rather different motivations, which seems reasonable, then we might assume 
independence and hence that 

Var[M F ] = L 2 (X B + X s )AtE[N 2 } . (16) 

2.4 The technical traders 

We consider now a second group of traders who only trade in response to a return 
created within the period under consideration. We shall model these trades in 
the same way as for the fundamental case, except now the trade arrival rates, 
instead of Ab,f are now /j, b (x) and Hs{x) where the buy and sell fii have the 
property that ^(0) = - we assume that in the absence of a price movement 
there are no technical trades. Treating the variation in trade size in an identical 
fashion, we have additional net buying pressure M T which is a random variable 
with 

E[M T ] - L{hb{x) - ^s{x))Atn , (17) 
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and a variance which again depends on the degree of correlation between the 
buy and sell actions. In this case we shall assume that the buy and sell actions 
are perfectly correlated, which makes sense for example if a trader is pursuing 
a mean-reverting strategy: 

Var[M T ] = L 2 (n B (x) - ^ s (x))AtE[N 2 } (18) 

It is natural to assume that the technical traders are operating independently 
from the fundamental traders and this will be done. 

2.5 The return impact function 

We introduce a log-return impact function I(q) that is a function with 1(0) = 
such that the aggregate buy and sell orders of both types create a log-return 
impact of the form 

Ax = X(M F + M T ) . (19) 

2.6 Summary of discrete model 

In the time interval At the return changes by 

Ax = T[M F + M T ) , (20) 

where T(q) is the return impact of a net order to buy q shares. The variable 
M F is a random variable with mean 

E[M F ] = L(\ B - X s )Atn (21) 

and variance 

Var[M F ] = L 2 (\ B + X s )AtE[N 2 ] . (22) 
The variable Mr is independent of M F , and has mean 

E[M T ] = L(fi B (x) - fx s (x))Atn (23) 

and variance 

Var[M T ] = L 2 (i2 B (x) - fi s (x))AtE[N 2 ] . (24) 

3 Linearization and the SDE 

In order to proceed to a continuum representation we shall now make some 
simplifying assumptions. These are 

1. Linearization of the return impact; 

2. Linearization of the \x functions; 

3. a Brownian motion model. 

We now consider each of these assumptions. 

3.1 Impact linearization 

For the return impact function all we know for sure is that no orders mean no 
price impact, i.e. 1(0) = 0. In general T(q) may be a very complicated function. 
We think it should not decrease as q increases. The structure of the order book 
may give it a staircase character 7 However, we shall assume that on certain 



Note we always work based on mid-prices so the central step is taken as absent. 
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scales we may linearize it - a staircase may look like a smooth line viewed from 
a "long way" away - in order to capture the gross feature of the model. So for 
some unknown constant uo we write 



The notion of a linear model of price impact was employed by Almgren and 
Chriss [3] as a component of both temporary and permanent price impact. It 
should be appreciated that this is very much in the same spirit as the assump- 
tions made in the Bakstein-Howison model [5], where the essential features of 
the order book are reduced to a spread and liquidity parameter pair, where 
liquidity is the reciprocal of the slope of the order book. This idea is a pow- 
erful one which may also be used to analyze derivatives, as considered e.g. by 
Mitton [30] . Thus liquidity is an implicit part of the model via the linearized 
price-impact. 

3.2 Technical trade linearization 

The appearance of trades depending on price movements within the time period 
under consideration will also have a complicated dependence. The existence of 
limit orders at various price thresholds will also create a staircase effect. We shall 
linearize this in the same way as the price impact, and make the replacement: 



where /x is an effective constant parameter that captures the gross slope of the 
function 8 We put in a minus sign due to the nature of limit orders coming into 
play to act against the direction of price movement. One would expect the effect 
of fj, to be negative on balance due to the effects of profit taking as well, unless 
the system is being overwhelmed by momentum trades. The appearance of the 
latter effect will be analyzed in detail later in this paper. 

3.3 Process approximation 

The trade arrival model is at this stage purely of Poisson type with a variable 
trade size. There are a number of ways in which this might be managed, de- 
pending on the details and frequency of trade arrivals. If there are a sparse 
number of large trades, this will be best viewed as a pure jump process. We 
might also have a situation where a large number of trades of moderate size 
together with sporadic large ones. This will generate in effect a jump diffusion. 
In the following simplification we assume that there are a large number of trades 
of moderate size so that we do not consider the jumps. The other cases will 
be investigated elsewhere. It is also possible, indeed likely, that the net price 
impacts of the two types of trade and their associated volatilities can be time 
dependent - this model docs not rule out stochastic volatility at all. However, 
in what follows we shall confine attention to constant parameters and the pure- 
diffusion view - even this subset of possibilities will demonstrate a rich structure 
through the emergence of a hybrid arithmetic-geometric stochastic process. 

In this context we now approximate the Poisson-type trade(s) arrival model(s) 
by independent Brownian motions centred on the mean arrival rate. We are lead 
finally to a discrete-time stochastic evolution equation for the return: 



Ax = uj x (M F + M T ) . 



(25) 



Hb{x) - ns{x) — > -fix x 



(26) 



Ax = ujL n[(X B - A s ) - fix]At + s x AWi + s 2 xAW 2 



(27) 



8 In both this case and the treatment of price impact we are not making differentiability 
assumptions and using a power series - the idea in both cases is that a complex staircase may 
be grossly idealized as a sloping plane. 
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where 

Sl = ^(X B + \ S )E[N 2 } , s 2 = y/fiE[N2] . (28) 
We now take the continuum limit and finally arrive at the SDE 

dX t = ( Ml - p 2 X t )dt + (TidWu + o 2 X t dW 2t , (29) 

where 

Pi = an(X B - A 5 ) , 
p 2 = apn , 

ui = ay /(X B + X S )E[N 2 } , (30) 
a 2 = a^pE[N 2 } , 

and a = Luj is the return impact of trading one lot of shares. The SDE given 
above is the basic description where we show the explicit contribution separately 
of the fundamental and technical trades. We can of course reduce it to an SDE 
with a single noise term as follows. If p is the correlation between the two 
Brownian motions, then we can write the SDE as 

dX t = {pi - p 2 X t )dt + ^a 2 + X 2 a 2 + 2paiX t a 2 dW t . (31) 

This is one of the class of "Pearson diffusions" considered by Forman and 
S0rensen [20]. The first detailed application to financial modelling that the 
author is aware of is the work by Nagahara [31]. It has a notable special case 
that we now consider. 



4 The Student equilibrium model 

A particular case of interest is obtained by considering pi = = p, so that we 
obtain the SDE 

2^2 



dX t = -p 2 X t dt +\J(Ti + Xlo\ dW t . (32) 

In the equilibrium situation, the quantile ODE [39] associated with this SDE 
reduces to 

d 2 Q f8Q\ ~ 2 2(a 2 + p 2 )Q 



3u 2 \ du) a 2 + a 2 Q 2 ' (33) 

Bearing in mind the results of [39] we see that we have a quantile function for 
a Student distribution with 

Q - r ^-— w{u) , (34) 
Verl + 2^2 

where w(u) is the standard Student quantile with degrees of freedom 

v= 1 + 2^| . (35) 

°2 

So it is clear that we need p 2 > for this to be a Student distribution. This 
of course corresponds to the requirement that the underlying SDE mean-revert 
to the origin, and this mean-reversion condition in turn allows an equilibrium 
to establish. This equilibrium origination of the standard Student distribution 
arises naturally in plasma physics [38]. The faster the mean-reversion rate is 
compared to the multiplicative volatility, the closer the system is to the normally 
distributed limit. The Student distribution also arises naturally in the modelling 
of asset returns [19, 36, 41]. 

In more generality we obtain a dynamic hyperbolic generalization of an OU 
process. This has been argued by Forman and S0rensen to be, in equilibrium, 
a natural candidate for a skew-Student model. In complete generality a still 
richer class of diffusions with or without jumps is possible. 
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5 The "hyperbolic O-U" SDE 

While we do not have a characterization of the full time-dependent aspects, 
some insight can be gained by reducing one form of the hybrid SDE to standard 
form. First we do some scalings to standardize Eqn.(31). We let 

X t = -^=Y t , S = o 2 JD . (36) 

Then the SDE for Y t is 



dY t = -^(l--)Y t dt+Z \ l+—dW t , (37) 



\ v J V v 

and this is a simple two-parameter form of the problem. We can re-cast this by 
setting 

Y t = ^sinh(Z t ) , (38) 

and this leads to 

dZ t = — tanh(Z t )dt + -jLdW t , (39) 
2 \ v 



or equivalcntly 



2 



dZ t = - ^-utanh(Z t )dt + cr 2 dW t . (40) 
With a time-scaling we have the non-dimensional form: 

dZ T = - ^vt<mh(Z t )dT + dW T . (41) 

The original variable X t is then given in terms of Z t as simply: 

X t =— sinh(Z t ) . (42) 

One can explore the full Fokker-Planck equation based on cither of equations 
(36) or (39). Equation (39) reveals the essential nature of the process: for 
small Z t and small times, the process is essentially OU in character. But the 
mean-reversion in the tails levels off and becomes much weaker. 



6 Pearson type IV: skew-Student equilibria 

We now turn to the more general case, where the SDE is 

dX t - (mi - lnX t )dt + ^al + X?ai + 2pa 1 X t a 2 2 dW t . (43) 

This time, from the results of [39] we obtain the equilibrium quantile ODE as 

d 2 Q (dQ\- 2 = 2[(pa l( T 2 - Ml ) + (o\ + M2 )Q)] 

dv?\du) (c7 2 1 +ajQ 2 + 2pa 1 a 2 Q) ' 1 ' 

This is then related to the logarithmic derivative of the density function f(x) 
as 

1 df(Q) _ 2\{po x o 2 - p^ + {a% + p 2 )Q)] 



f(Q) dQ ~ (a 2 + a 2 Q 2 + 2pa 1 a 2 Q) 
We can solve this ODE and find that, after careful normalization, 

2 -, -(f+l)/2 



fix) = k 



1 + 



x — A 
a 



cxp 



-v 2 tan 1 



x — A 
a 



(45) 



(46) 
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where the parameters are given by 



CT 2 



0-1 

CT 2 
M2 



m = 1 + 



1 + 2 



M2 
^■2 ' 



(47) 



2(/iio- 2 + pcriM2) 



fc = 



en of 
r( 



2 , 



r(^) 



This is the class Pearson Type IV distribution, which is one candidate for a 
choice of "skew-Student" distribution, with a rich variety of skewness and kur- 
tosis in the structure. A useful guide to the properties of the Type IV Pearson 
is given by Heinrich [23], who uses the parameter m above, and his v is our v 2 , 
i.e. his density is 

2 



f{x) = k 



1 + 



A 



exp 



-z/tan 



A 



(48) 



Transferring the results of [23] to our own notation 9 we can identify the mean, 
provided v > 1, as 

E[X] ' ^ 



A 

v - 1 

The variance exists provided v > 2 and is then 

a 2 ((^ 



S 2 = £[X 2 ] - E[X} 2 = 



i) 2 + 4) 



(49) 



(50) 



iv-lYiv-2) 

The third moment can be calculated provided v > 3 and leads to the normalized 
skewness as 

E[iX-E[X]f] 



£ 3 



u-3 



iv - l) 2 + vl 



(51) 



The fourth moment exists provided v > 4 and may be expressed through excess 
kurtosis, which is 



E[(X - E[X]f 



-3 = 



6 (z/i - 3) (i/i - 1) 2 + 6 (5^i - 11) v\ 



(52) 



(^ 1 -4)(i, 1 -3)((^ 1 -l) 2 + ^) • 

When ^2 = the skewness is zero and the excess kurtosis reduces to the well- 
known expression for the pure Student distribution: 6/(i> — 4). We verified 
the translation of these expressions to our notation by the computation of the 
moments by simple numerical integration for numerous parameter values. 



7 Towards the full dynamics 

The achievement of an equilibrium is not realistic for most trading periods, 
especially during a panic of the credit-crunch period. The pure equilibrium 
analysis above is meant to indicate how a rich variety of distributional types may 
emerge from a simple model, and no more. We now turn to the full dynamics. 
We can consider this from both the SDE point of view and in terms of the 
Fokker-Planck equation for the time-dependent density. 

9 Users of the standard student 't' are perhaps more used to working with the degrees of 
freedom parameter v. 
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7.1 Formal solution of the SDE 

The solution of the SDE can be given explicitly, and gives a structure similar to 
the functional of pairs of Levy processes analysed by Carmona, Petit and Yor 
[11], in [6], where the structure arising from a pair of Levy processes is regarded 
as a generalized OU process. Here we will give the explicit integration of the SDE 
in the case of non-zero correlation for the Brownian motions via a correction to 
the drift somewhat similar to that employed for quanto options. The following 
version of the argument is based on the idea of treating the SDE in much the 
same way as one treats the OU SDE, by using an exponential integrating factor. 
In this case the integrating factor is itself a stochastic quantity. 10 Starting from 

dX t = (mi - (J,2X t )dt + (j-tdWu + o 2 X t dW 2t , 

we introduce the stochastic integrating factor 

I t = exp[-(7 2 W 2t + (/i 2 + iof)t] (53) 

and observe that 

dl t = (M2 + c\)hdt - <7 2 I t dW 2 t . (54) 
Now we apply the integrating factor by defining 

Qt = XJt (55) 

Then application of Ito's lemma and a short calculation gives 

dQt = - po\u 2 )dt + oidWtt) (56) 

and by integration, with X = = Qo, we obtain 

Qt = [ Isdm ~ P<Tia 2 )ds + (JidWis) (57) 
Jo 

X t = f (Ir 1 Ia){(Hl-(XTl<T2)d8 + (ndW le ) . (58) 

Jo 



and hence 



Now we observe that 



It 1 I s = cxp [<j 2 (W 2t ~W 2s )+ (/x 2 + l -al) (s - t)] (59) 

and so by a time reversal u = t — s, with associated time-reversed variables W i} 
we can write 

X t = J ((/ii - p<Ticr 2 )du + cridWiu) exp [<j 2 W 2u - (p 2 + -<J 2 )u] . (60) 

This reveals that the solution is the integral of one Brownian motion against 
the exponential of a second Brownian motion. Thus the process is a clear 
generalization both of the Asian SDE (cti = 0) and the OU SDE (cr 2 = 0), with 
the resulting complications. We can further reduce the expression in terms of 
the degrees of freedom variable v, to 



x t= ((Mi - p(Ti<J 2 )du + o x dW lv ) exp[a 2 W 2u - -cr 2 u] . 
Jo z 



(61) 



10 I am grateful to D. Crisan and M. Yor for pointing out the method for the zero correlation 
solution - personal communication, Kyoto 2009. 
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For future reference we note that in the case v = = p = p,\, we have 

X t = u x f dW lu exp[<r 2 W 2u ] . (62) 
Jo 

so that X t is just the integral of one pure Brownian motion against the expo- 
nential of another. By comparison with the hyperbolic OU SDE obtained in 
Section 5, we can observe the recovery of Bougerol's identity [8, 44], that in this 
case X t is distributionally the sinh of another Brownian motion. We shall see 
this special case emerge again when we solve the Fokker-Planck equation. 

7.2 Dynamic moment evolution 

A partial characterization of the full dynamical situation may be given in terms 
of the evolution of the moments. Let us set 

e n = E[X?] . (63) 

Then elementary analysis gives us the sequential families of ODEs: 

^ + (M2™- ^n(n-l)<r|)e n = ^n(n- l)cr?e„_ 2 + {p,\n + n{n- \)po\G 2 )e n -\ , 

(64) 

with eo = 1. The evolution of the mean e\ is therefore governed by 

-^+P-2e\=p,\ (65) 
and, as in an ordinary OU process, evolves according to 

d = X e-^ + ^(1 - er^ 1 ) . (66) 

M2 

With our conventions X n — at the start of the trading period so we then have 

ei--(l-e^ 2t ) , (67) 

M2 

which will settle down to ^1/^2 if ^2 > and grows exponentially otherwise. 

7.3 The explosion of variance 

Solution for the higher moments is straightforward but leads to rather unwieldy 
formulae in general. In the special case where p = = fi\ the variance V(X t ) 
may be written in the tractable form 



V(X t ) = 



e -al(v-2)t 



ob - \<Ao\{v - 2)t 2 + 0(t 3 ) , (68) 



where, as before v = 1 + 2^2/ 'o\. Once the market has kicked off (the behaviour 
near t = being always Gaussian), the market dynamics are thus critically 
dependent on the sign of v — 2. If the strength of the mean-reverting trades is 
such that p 2 > cl/^ the market settles down. If instead p 2 < cf/^ the variance 
grows exponentially. Note that there is a region < p, 2 < a 2 /2 where the 
average level stabilizes but the variance does not. If [i 2 < both the average 
level and variance grow exponentially. If v > 2 the distribution settles down 
to the Student equilibrium already analyzed, the with the condition v > 2 
guaranteeing a finite variance. 

The exponential growth in variance is not a new idea - it has been present 
in the price evolution model of geometric Brownian motion for decades. What 
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is being suggested here is that technical trade effects in the evolution of the 
Zog-return lead to a variance explosion in the log-returns due to these being a 
hybrid arithmetic-geometric hybrid. 

We can now return to the famous 25a issue. If one's perception is that 
normal Gaussian behaviour is to be expected, with the classical variance a\t 
in the log-returns, then the inclusion of technical trading effects via the hybrid 
model causes the actual variance to differ by a ratio, that we call the variance 
explosion factor 

V(t) 



v E (t) = 



a\t 



(69) 



In general the variance explosion factor may be found by solving the ODE for 
e 2 . In the special case considered above we have 



v E (t) 



1 



a\t{v - 2) 



1 - e 



(70) 



If the markets arc minded to settle to equilibrium this ratio tends to zero. But 
if v < 2 we have instead 



V E (t) = 



1 

at 



e (°.t) _ i 



<jI{2-v) > 



To return to the likelihood of 25cr events, we note that if 
a > 6.4746 , i.e. a\{2 - v) > 6.4746 , 



(71) 



(72) 



then V E {t) > 100 and a 25c event based on an initial perception of variance 
aft is no less likely than a 2.5ct event with the right variance (and with a 
different, non-Gaussian, distribution). Such events may then occur repeatedly 
if the technical market effects are strong enough, even without incorporating 
the additional effect of net price pressure due to fundamental trades (hi ^ 0). 

We see that a variety of different dynamics are possible, with quite small 
shifts in the mean-reversion strength of technical trades making a dramatic 
difference to the return distribution. 



7.4 Dynamic distributional aspects 

Having identified the impact of the arithmetic-geometric hybrid on the moments 
we now need to understand the full shape of the distribution. If the distribution 
is significantly fat tailed then this can also amplify the likelihood of extreme 
movements. There are three ways of proceeding. One approach is to transform 
the problem to the form analyzed by Wong [43]. The second is to seek a fresh 
approach to the Fokker-Planck equation. The third is to exploit the represen- 
tation of Section 7.1 to obtain further insight from a conditionally Gaussian 
representation. The representation given by Wong involves the integrals of hy- 
pergeometric functions, so further insight is helpful, both to understand the 
theoretical behaviour and to provide more practical tools for practitioners. 

First we should note that the problem of finding the time-dependent PDF 
in the case p — = /xi can be written down by transforming the problem to 
that analyzed by Wong [43]. This is the form also used by Nagahara [31]. In 
our notation, Wong's solution has p = = p,i, v = 2a, a\ = a 2 = \pl. The 
Wong paper from 1963 is available on-line from E. Wong's web page at 

http : //www . eecs . berkeley . edu/~wong/wong_pubs/wonglO . pdf 

An issue for financial analysis is the detailed practical implementation of 
this density and the extraction of insight. Wong also notes that the solution 
becomes much more straightforward in the case when a is a positive integer, and 
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gives Gaussian functionals for those cases. In general however, the density is an 
integral of a hypergeometric function of complex parameters. Such expressions 
are tractable in an advanced computing environment such as Mathematica, but 
present challenges in C or other less mathematical languages, in much the same 
way as for Asian options. Our approach here will be to obtain some insight 
beyond the work of Nagahara and Wong by direct analysis of the Fokkcr-Planck 
equation and relating the outcomes to the theory of integrals of exponentials of 
Brownian motion, and to the Legendre equation common in physics. We will 
see that matters are in fact more straightforward when a = u/2 is any integer, 
and this is a manifestation of the discrete symmetry of the Legendre equation. 
Negative cases correspond to momentum-dominated markets, and positive cases 
are mean-reverting and asymptotically Student. 

A further point of practical implementation of the Wong formula is that 
given in general it involves complex integrals of hypergeometric functions, it is 
not necessarily of greater ease of use than inversion of the transform that we shall 
derive here, which is also hypergeometric. The underlying issue is whether it is 
helpful to separate out the discrete eigenvalues of the Sturm-Liouville operator 
from the continuous ones. The Wong formula has this separation - the transform 
representation here does not. Both require an inversion by complex integration 
in general. The derivation from first principles of the transform given here al 
least offers an alternative, and indeed we can transform it into different forms to 
make the analysis of positive and negative v more transparent. It is hoped that 
the options for representation here will stimulate further analysis. The matter is 
somewhat similar to the older analytical studies of Asian options triggered by the 
Gcman-Yor model [21]. Direct inversion of the transform in Mathematica was 
possible with a few lines of code [35], as long as the parameter a 2 T was not too 
small. The work by Linetsky [25] established the useful spectral representation, 
but probably the simplest solution to the low volatility evaluation problem is 
an asymptotic analysis [13, 46]. 

To get at the PDF in general we must analyze the full Fokker-Planck equation 
in the form 



df(x,t) 
dt 



d_ 

dx 



1 d 



-(fii-fi 2 x)f(x,t) + - — [(al+x 2 al+2paia 2 x)f(x,t))] 



(73) 



with the initial condition f(x,0) = S(x). One route to this is via the Laplace 
transform with respect to time. So let 



f{x,p) 



f(x,t)e- pt dt 



(74) 



Then the Fokkcr-Planck equation gives us, suppressing the independent vari- 
ables. 



pf-S(x,0) 



dx 



1 d 



(fix - fjb2x)f + - + x 2 a\ + 2pa 1 a 2 x)f)] 



, (75) 



This is now a Green's function computation on the transform. For x > and 
x < we need two independent solutions of 



Pf 



d_ 
dx 



Id 

-(Mi - H2x)f +2 0% + a;2<J 2 + 2p<Tif7 2 a;)/)] 



(76) 



with the junction condition that / is continuous at x = 0, and a jump in the 
first derivative. This condition, integrating about zero, is 



<J/ (o+,p)-^(o-,p) 



dx 



T 2 ' 



(77) 
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7.4.1 The Gaussian case 

In order to make this approach clearer, we first pursue the standard Gaussian 
problem where <j 2 — — /i 2 - The solution of the transformed ODE vanishing 
as x — > ±00 and with the correct junction condition at zero is 



e fj,ix/af (expr-^x/u? + 2poV\ if x > 0, 
1 ' y/fi + \exp[+^ VM? + 2p<7?] if x < 0. 
Inversion of the two cases leads to the single well-known formula 

f(x, t) = — L= expf-Or - ^t) 2 /(2a 2 t)] . (79) 

\J2-K(j{t 

7.4.2 A dynamic Student distribution 

The next case of interest is when p = = for which we previously demon- 
strated a Student t equilibrium under certain circumstances. The Laplace trans- 
form of the Fokker-Planck equation is, for x ^ 0, 

l -(o\ + a 2 x 2 )f"(x,p) + ( M 2 + 2a 2 2 )xf'(x,p) + ( M2 + o\ - p)f(x,p) = . (80) 
This equation may be simplified somewhat by setting 

f(x,p) = (a 2 + a 2 x 2 )^ 1+ ^g(x,p) , (81) 
and the ODE for g(x,p) is then 

(a 2 + o\x 2 )g"{x,p) - 2x f i 2 g , (x,p) - 2pg(x,p) = . (82) 

We have already worked out the equilibrium case when p = and g is constant 
in x. The management of such an equation is straightforward in the special 
case \i 2 = —(T2/2, as discussed in [33]. We shall use the change of independent 
variables indicated in [33] to treat the general case, and indeed this is almost 
the same change of variables that took us to the hyperbolic OU picture. We 
introduce z(x) with the condition that 



dz _ 1 

dx ^/a 2 + a\x 2 

and fix the arbitrary constants so that 



(83) 



z=-sinh- 1 ^ N ) . (84) 



C 2 V 01 

Our equation for g expressed in terms of z is then just 

%\ - (2fx 2 + of) — tanh( ( r 2 z) ^j- — 2pg = , (85) 
az a 2 az 

or in terms of the degrees of freedom parameter v = 1 + 2fj, 2 /a 2 , 

— va 2 tanh(cr 2 z) ^~ - 2p.g = • (86) 

We previously remarked on the Bougcrol identity [8] for the case when v = 
0. This finds a simple expression here (see also Section 3 of Matsumoto and 
Yor[28]). In this case the transform is 

f( \ - 1 1 JVv^'W if z, x > 0, 

" + ^ \e+v*P*(*) i£z,x<0. ( ' 
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and the inversion gives us 



f(x,t) 



1 



: exp 



1 



2o\t 



[sinh 1 (<7 2 x/(T 1 )Y 



^/2-Kt{o\ + o\x 2 ) 

which is the density arising from a change of variables z — > x on the z-density 

1 



exp(-z7(2t)) , 



(89) 



and is a significantly fatter-tailed object (while in the neighbourhood of x = re- 
sembling a Gaussian with variance cr 2 t\) The presence of significant momentum 
trading has fattened the tails in this case. We shall now explore the properties 
of the PDF given for the special case hybrid given by Equation (76). 

7.5 The hidden menace 



6 - 



4 - 



.2 - 



Figure 1: PDFs for Gaussian and special case hybrid, t = 0.1. 




Figure 2: PDFs for Gaussian and special case hybrid, t = 1. 



Even the very special case solved for exhibits some interesting features. At 
the start of the trading period the distribution of log-returns is barely distin- 
guishable from Gaussian, as shown in Figure 1 for the parameters o\ = o-i = 
l,t — 0.1, where the hybrid is overlaid with the Gaussian. 
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Figure 3: PDFs for Gaussian and special case hybrid, t — 5. 



However, as time passes the hybrid distribution spreads out more, in a man- 
ner consistent with the variance explosion formula. In Figure 2 and Figure 3 
we show the hybrid PDF overlaid with the Gaussian at times t = 1 and t = 5 
respectively. At all times the probability of a very small move remains at the 
Gaussian level - the PDF osculates the Gaussian at the origin. The overall be- 
haviour represents the hidden menace of these processes. It starts off looking 
Gaussian with variance aft; the probability of a very small movement remains 
near the Gaussian value, yet dependent on the size of a 2 the probability of 
extreme movements grows exponentially in time. 



8 More general solutions 

We now look at the case where v > 0. We go back to Eqns. (84-86), and change 
independent variable to 



u = sinh 1 (o'2^/ci) = 022 



(90) 



Setting s = 2p/cr|, and letting ' now denote d/du, Eqn. (86) may be reorganized 
as 

(e u + e- u )(g"(u) - sg{u)) = v{e u - e- u )g\u) . (91) 

Previously we considered the case v = 0, where we could write down a solution 
decaying as x — » +00 as a single exponential in e~^ u . When v ^ we have to 
proceed differently, but the presence of the two exponentials in the coefficients 
suggests a solution approach. To tidy up a little we make a further change of 
variable to 



So with u 

as 





w 


— log 


uu we 


( l - 













-wd/dw and our ODE may be rewritten 



d 2 9 



dg 



dw 2 dw 
We seek a power series solution in the form 



dg 

1 — 
dw 



(93) 



E 

fc=0 



a k w 



fe+7 



(94) 
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After some standard manipulations we find an indicial equation in the form 

7 2 + - s = , (95) 

and the required root, for v > 0, to get the right behaviour as x — > oo, w — > 0, 
is 

^=\A+T"^ (96) 

The resulting recurrence relation simplifies to 

a k+2 {k + 2){k + 2 + 2 1 + v) = -a k {k-v){k + 2 1 ) . (97) 

After some experimentation with hypcrgeometric functions and some work with 
Mathematica we are able to recognize the solution in the form 

g = a {p)w-< 2 F! (7,-^7+ ^ + 1;-™ 2 ) . (98) 

This representation of the solution has the nice property that we can see that 
the hypcrgeometric function reduces to a polynomial if v is an even integer. The 
case v = has already been exhibited. Before discussing other such simple cases 
we must complete the solution and determine do- This involves the application 
of the jump condition on the derivative at the origin, assuming an even solution. 
After some algebra we find that 

ao(p) - Z 1 ; , (99) 
cr 2 S2(z/,7) 



where 



(100) 



After some use of Kummer's identity and variations (specifically identities 15.1.21 
and 15.1.22 from [1]), we are lead to 

r (2) r (2(7 + ^ + 1)) 
We arrive at a closed form for the Laplace transform of the density as 

- = a^nPT (?) r (|( 7 + v + 1)) 2*1 (7, -f ; 7 + | + 1; -^ 2 ) 
A ' w ^a 2 r (7+I + 1) K + x 2 °D 

We remind the reader that in the use of this expression, 

. 2 P 



(102) 



w 











I/ 2 


V 


\l s + 


T 


~ 2 ' 


( °2X 

V CT i 


+ 1 





(103) 



(7? 



It is now easy to check the known special case, when v = 0, for the hypergeo- 
metric function simplifies leading to 
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We now also have a new family of relatively simple cases when v is an even 
integer. For example, when i/ = 2we have 

f(x,p) = - G } (— + , 7 = ^^TT - 1 . (105) 

2cr 2 (erf + 2H0-2) V 7 7 + 2/ 
In the case v = 4we have 

'(3 + 7)w 7 2 W t+ 2 ( 7 + l)wT+ 4 



f{X,P) Aa 2 (a 2 + W 2 ) V 7(2 + 7) ' 7 + 2 ' ( 7 + 2)(7 + 4); ' (106) 
7= Vs + 4 - 2 . 

8.1 A Legendre representation 

The solution for the Laplace transform when ji\ = p = has now been charac- 
terized. However, the representation above is not the only one that can be given. 
The solution is in the form of a hypergeometric function 2-F\(a, b, a — 6+1, — w 2 ), 
where a = 7, b = —v/2 in this case. Such a pattern of arguments to 2-F1 leads 
to several other equivalent representations. For example, five others are given 
in [1]. One of these is of particular note, and is associated with identity 15.4.15 
of [1]. This is the relation, valid for —00 < z < 0, 

2 F l {a,b,a-b+l,z) = T{a-b+l){\- z)-\-z) b l 2 - a l 2 Pt- b a (^^j , (107) 

where P^f is the associated Legendre function. Making the relevant substitu- 
tions with z = —w 2 and simplifying the result gives us another formula for the 
transform: 

/>,„) = ^^r(|)r(l±^)(cosh U )+^+ 1 )p; / f ^(Itanh.l) , 

(108) 

where u = sinh 1 (ct 2 x/cti). This gives us a nice interpretation of the result. 
The quantity L = v/2 is the quantity in physics normally associated with the 
angular momentum of a quantum system such as the hydrogen atom. The 
Legendre functions in that case arise naturally via separation of variables of the 
Laplacian operator. The case L an integer is particularly simple as then the 
associated Legendre functions can be written in terms of polynomials. While 
these are not the ordinary Legendre polynomials (the quantity M here is a 
complex transform variable and not an integer) they are well known - see e.g. 
[48] for explicit forms for the first ten. Armed with the hindsight of Eqn. (108) 
one can now revisit the entire solution process and make a change of variables 
to reduce the transformed Fokker-Planck equation to Legendre's equation. By 
making the change of variables 

the transformed Fokker-Planck equation can be written as 



v,v v _ 2p/a 2 + u 2 /4 

2 [ 2 + 1 (l-g 2 ) 



%) = (110) 



which is the Legendre equation with parameters L — v/2 and M 2 — 2p/a 2 + 
i^ 2 /4. Such equivalent expressions reveal to us that the transformed distri- 
bution is a well-known mathematical object, and perhaps allow evaluation in 
computation systems where a full implementation of 2-F1 might not be available. 
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However, for further analysis here we will work with the original hyper geomet- 
ric system and the Legendre representation. The hypergeometric form is well 
adapted to an analysis of the tails, due to the power behaviour in w, and there 
are hints that this variable may be more useful in the full general case where 
fi± =/= 0. There are a limited number of known inversions of Legendre functions 
where the transform variable is an index to the function and these are also under 
investigation. For now we turn to another special case for more detailed study. 



8.2 A "chameleon" distribution 



In this sub-section we shall introduce a one-parameter family of dynamic distri- 
butions with the following interesting properties: 

• The distributions arise from stochastic differential equations; 

• The mean is identically zero; 

• The variance is of the form o\t\ 

• The behaviour is initially Gaussian, in standard form. 

• The distribution tends to a non-Gaussian steady-state. 

This is just a matter of specializing the analysis above to the case v = 2. Writing 
out the solution for the transform more explicitly, we have 



f(x,p) = 



2<t 2 (<j 2 + ct 2 x 2 ) 3 / 2 



-[v'i+T-llKx) 

VsTT- l 



+ 



-[V^+l+l]|«(x)|] 



where 



2p 



, u(x) = sinh 1 (ct2x/cti) . 



(Ill) 
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This may be inverted in closed form, making careful use of identity 29.3.88 
from [1] and some standard Laplace transform identities. After some careful 
simplifications we are lead to the following density function: 



+ 



r u(x) 2 

<j\ exp[ — - 



2 



CT 2 cr 2 



2 (a 2 + x 2 a 2 ) 3 / 2 



\u{x)\+ta\ 



0-2 



\u(x)\ — t&2 
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where $ is the standard normal CDF. This is the probability density function for 
our "chameleon distribution" . We have obtained it by solving the Fokker-Planck 
equation from an SDE. Its mean is zero and is variance satisfies 



V(X) 



ah 



Vct 2 , t 



(114) 



Its asymptotic behaviour as t — > is obtained by considering just the first 
exponential part of the expression, which we see tends to 



r u(x) 2 

a i exp[ 



2t<7 2 



exp[ 



Itai 



fWt (<r 2 + x 2 al) 



(115) 



where the last approximation arises as the support of the distribution contracts 
about the origin, allowing us to expand the arcsinh and denominator. So the 
distribution starts off in standard Gaussian form. For i^oowe just note that 
the first term tends to zero and the second line tends to 
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which is the density of a scaled Student t distribution with two degrees of 
freedom. This of course has infinite variance. We have therefore demonstrated 
the list of conditions claimed at the start of this Section. Parametrized by 02, 
there are infinitely many SDEs of the form 

2 , 

dX t = —±X t dt +\]<ri + ^XfdWt , (117) 

whose variance is given by the standard formula aft, normally associated with 
the simple case 

dX t = a x dW t ■ (118) 

This is a useful reminder that linear evolution of variance does not in any way 
imply the elementary Brownian model. Although we have produced a rather 
exotic density function, the underlying dynamics as given by Eqn. (117) are 
very simple. 



8.3 Momentum-dominated markets 

We now consider the case v < 0. We can go back to the Fokker-Planck equa- 
tion (80) and seek more appropriate changes of variables, or explore the con- 
tinuation of the v > procedure to v < 0, and observe that it must remain 
valid. Equivalcntly we can exploit the symmetry of the Legendre representation 
under L — > — (L+ 1) to derive the appropriate transform. Perhaps the simplest 
approach is to exploit a hypcrgcometric identity (see e.g. [1]), that is also an 
encoding of the Legendre symmetry: 

2 F 1 (c -a,c-b,c,z) = (l- zf a+b - c \F 1 {a 1 b, c, z) (119) 

We apply this identity with a — 7, b = —v/2, c = 7 + v j2 + 1, z = —w 2 

2F 1 (" / ^, 1 +^ + l,-w 2 ) = (l + w 2 )^ +1 \F 1 (^ + l r/ + is+l, 1 +^ + l^w 2 ) 

(120) 

Next, given the form of w, we note that 



\ + u , 2 = 2 ^^a 2 + a 2 x 2 (121) 



and putting these observations together with the transform density representa- 
tion Eqn. (102) we obtain a representation better adapted to v < 0: 

- _ 2 7+^ + i r (2) r (1(7 +.+ !)) 

V^CTiCT 2 r (7+ § + 1) ( 122 ) 
X 2^l(^ + l,7+^+ 1,7+ ^ + l,-w 2 ) . 
where we also recall that 



7= V S+ 4 2 

We now see that another set of special and simple cases emerge. For v — — 2 
we have, simplifying, 

f(x,p) = — - = — . (123) 

0102(7 "I J 010-2VS+1 

For v — —4 we have, simplifying, 

J ~ {X ' P) = 2a~^ (7^1 + ~i) (124) 
7 = Vs + 4 + 2 
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And for v = — 6 we have, simplifying, 

- 1 / ( 7 -2)w^- 5 2«;T- 3 ( 7 -4)wT- 1 

4o^l v ( 7 -5)( 7 -3) ( 7 -3)( 7 -l); (125) 

7 = Vs + 9 + 3 

These of course mirror the three cases v = 0,2,4 described earlier and exhibit 
the Legendre symmetry, but without the Student denominators, as is natural 
for the increasingly explosive market behaviour that these new cases describe. 

8.4 Generalized Bougerol Identity 

The case v = —2 merits further study as it is the simplest momentum-dominated 
case, and it turns out links to some nice insights on hyperbolic Brownian motion 
[2], as also summarized in [9]. The x-space density is obtained by inversion as 



'., , exp( 1' '1 > (126) 



where 

• 1,-1 U2X 
u = smh 

Now transforming to the w-density, with dx — > ^ cosh(u)du, we obtain the 
u-space density as: 



9 {u,t) = ^CO S h(u)f( X ,t) = _A^( e Hu-«ltf/{2«tt) +e - { u + altf, { 2*lt)\ 
02 (T2\/27rt V / 

(127) 

so that in the dimensionless time coordinate r = u\t we see that the variable 
<JiXtj(j\ is distributed as 

sinh(T^ T + Yt) (128) 

where Y is Bernoulli and W T is a standard Brownian motion. This corresponds 
to exponential drift +1 when viewed as the exponential of a Brownian motion. 
It is possible that similar simple interpretations exists for drift +2 (y = —4) 
and indeed for the simple cases considered earlier v = 2, 4 where the densities 
also reveal a similar superposition of opposite and multiple drifts. 

Perhaps the most significant practical consequence is the emergence of a 
bi-modal density, where there are essentially two densities moving apart at a 
constant speed. In the absence of fundamental drift this makes good sense for a 
momentum-dominated model, and is expected to be a generic feature for v < 0. 

8.5 Conditionally Gaussian form for p = 

The third route to characterizing the density involves the formal solution of the 
SDE in the form 

X t = / (^du + a!dW lu ) exp[a 2 W2u - ^crju] . (129) 
Jo z 

This is valid for zero correlation between the two Brownian motions. If we 
condition on the exponential Brownian motion 11 the conditional distribution of 
X t is readily seen to be Gaussian with mean m and variance v, where 



m - 



m J ducxp[a 2 W2u - 2 CT 2 M ] (!30) 



: I am grateful to D Crisan for this idea. 
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v = o\ I exp[2<72W2u - vu\v\ . (131) 
Jo 

This is an elegant representation that makes it clear that the mean and variance 
are Asian- like quantities. This also opens up the use of other exact and approx- 
imate distributional representations - this will be pursued elsewhere. However, 
we note that this is a natural outcome if we wish to have a conditional Gaus- 
sian model whose total distribution is asymptotically a Student t, in that the 
conditional variance must asymptotically be inverse Gamma. This condition 
is satisfied by a conditional variance given by an Asian-like variable, i.e. the 
integral of the exponential of a Brownian motion. That this is created by a 
simple price-feedback model we think is interesting. One final point is that this 
type of representation is frequently used in the probabilistic literature. For ex- 
ample, the description of the Bougerol identity as phrased by Yor [44] is in the 
equivalent representation of a Brownian motion with an Asian time change, e.g. 

sinhB t ( L t^) 7At (132) 

where B and 7 are independent Brownian motions and A t is an Asian form. 

9 States of the market 

Having looked at some typical dynamics, we look at the overall picture. Our 
model has four parameters, and we now give them names: 

• o 1 , the fundamental volatility; 

• 02, the technical volatility; 

• hi, the fundamental drift; 

• fi2, the technical drift; 

While we have not incorporated fii m to any detailed analysis thus far, its in- 
terpretation is clear. The main influences on the state of the market are the 
parameters fi2, 01,0-2. In fact, it is the balance between these parameters that 
matters. A critical quantity is 

i/=l+^J. (133) 

If an equilibrium is achieved, this is the degrees of freedom of the associated 
Student distribution that results. But now we see that it plays an essential 
dynamical role: 

• if v < 2 the variance explodes exponentially; 

• if v = 2 the variance remains in Gaussian form, but any member of the 
chameleon family may exist; 

• v > 2 the variance tends to a constant. 

The circumstances under which the distribution attains an equilibrium are sub- 
tly different. We know that when v = 0, m = 1/2 the solution is eternally 
dynamic and the PDF has been calculated explicitly. When m > 1/2, > 0, 
the equilibrium solution exists and has a normalizable PDF in Pearson IV form 
[23] . So there is a range < v < 2 where the equilibrium exists but the variance 
explodes. It would be useful to get a better grip on the dynamic Cauchy case 
that sites in the middle of this zone with v = 1 = m. This case corresponds 
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to fj,2 = 0. When ^ < 0, and the market is momentum-dominated, matters are 
completely different and the distribution may become bimodal, as exemplified 
by the Alili-Dufresne-Yor generalization of the Bougerol identity for v = —2 [2]. 

The volatilities themselves play a different role. Inspection of the hyberbolic 
OU equation indicates that it is a 2 that sets the time- scales, and the ratio 
(j 2 / (Ti determines the scale on which asset price movements are affected by the 
price feedback. So in the absence of fundamental drift the market state is best 
characterized by the triple: 

• v = 1 + ^£ - determines the market condition; 

• <j 2 defines the time-scale; 

• 02/01 defines the asset price scale. 

Since all but the first are scaling variables, we can observe that when \i\ = 0, 
the parameter v is critical. Historically this has been the elementary 'degrees 
of freedom" parameter for static Student distributions. It now plays a critical 
dynamical role. If we set r = <j\t the underlying dimensionless SDE is 

dZ T = - V - tanh(Z T )dr + dW T . (134) 

with all other parameters removed by transformation. We have explicit time- 
domain solutions for the resulting PDF for v = 0,2 and the Laplace transform 
for other values. 



10 Initial VaR analysis 

The role and controversy around VaR (Value at Risk) have been discussed ex- 
tensively elsewhere. There is a fundamental issue of principle about estimating 
the probability of extreme events. We never know when there is a "tsunami" - 
type event that falls outside all of our traditional modelling. However, within 
any proposed model we should also look at estimating the risk numbers, in order 
to at least give an improvement based on the analysis we have been able to do. 



10.1 Simple Gaussian VaR 

Consider a percentile point u in the unit interval < u n < 1. Let Q denote 
the standard Gaussian quantile function associated with the CDF <f>(x). So we 
have the identity 

*[Q(«o)] - u . (135) 

In a simple analysis we might take u = 0.25, for example. In a standard 
Gaussian model, in the absence of fundamental drift, the signed Var would just 
be 

sVcltq = aiVtQ(uo) (136) 
and the VaR would be the absolute value of this scaled Gaussian quantile. 



10.2 Simple hyperbolic VaR via Bougerol 



The simplest more general case is that we have previously considered, which 
is the v = market state. Bougerol's identity gives us the density, and the 
cumulative distribution function is easily computed from the density function 
in Section 7.3.2, as 

1 . , _1 /<J2V\ 



F(y) = $ 



CT 2 



sinh _1 (- 



0i 



(137) 
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The inversion of this relation for the quantilc is now trivial, and we find that 



The value at risk is then seen to be hyperbolically greater than the Gaussian 
case, and the only coincide in the absence of any technical volatility or for very 
small times. 

10.3 More general hyperbolic VaR 

It is an interesting question to pursue the VaR estimates within the most general 
model. The Bougerol identity and its speculative generalization might suggest 
that there is an approximate approach for the momentum-dominated case by 
adjusting the simple hyperbolic VaR by a drift correction. In the mean-reverting 
case one might consider using the stationary asymptotic form. This is particu- 
larly simple for the Student case, where there is a lot of developed technology 
for the T-quantile. There are closed forms for v = 1,2,4 - see [36], and the 
on-line entry at [47], and series solutions for positive real v to any required 
precision [39]. Bailey's elegant generalization of the Box-Muller formula and 
its polar form [4] allows rapid adaptation of many embedded risk management 
models. A systematic mathematical treatment of an optimal VaR formula for 
the general case will be pursued elsewhere. 

11 Conclusions and further work 

We have developed a simple hybrid trading model that mixes both fundamental 
and technical trades. The general form of the model allows for jumps and 
Brownian motion together, but the theory here has been more fully developed 
for the pure Brownian case, where we have obtained several insights into the 
model proposed by Nagahara [31]. We obtain a hybrid process for the log-return 
that is a composite of arithmetical and geometric types, and that is capable of 
exhibiting a variety of behaviours depending on the relative strength of the 
fundamental and technical components. 

Even when the market settles down we obtain a rich family of distributions 
including fat-tailed Student and skew-Student models as a special case. It is 
therefore possible to attribute some of the skewness and fat-tailed behaviour in 
asset returns to a simple composite model of traders acting on different logic. 

In the fully dynamic context a the likelihood of extreme events is greatly 
increased by the use of a hybrid arithmetical-geometric process for the log- 
returns. Such hybrids can exhibit the phenomenon of variance explosion at the 
same time as hiding their non-Gaussian nature very effectively. Work on the 
full time-dependent form, with all parameters non-zero, is in progress, as is the 
VaR and related analysis. 

Another thread of research involves the use of multivariate extensions of 
these models, where the coupled distributions arise through coupled diffusions 
of Pearson type. In this way multivariate distributions with marginals drawn 
from the Pearson family may be constructed, as discussed by Shaw and Munir 



The adaptation of the model here to full risk-neutral pricing requires further 
work. The production of a density for the case where the drift is constant needs 
to be undertaken. But there are two further issues. First, the model here is 
regarded as adapted to short-term return modelling and the risk consequences. 
In reality different trading periods will have different market states. Some days 
will settle due to dominant mean-reversion, while others will have explosive 
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behaviour. So it is not regarded as appropriate to simply replace the log-normal 
distribution by the ones implied here - a more subtle approach is needed, though 
one could consider short-term options within this model. A further difficulty 
is whether the dynamic hedging process remains viable in the presence of the 
potentially explosive movements within this model in a momentum-dominated 
state. While it is mathematically possible to define the hedging process, in 
reality the explosion of variance could result in rapid movements that render it 
infeasible. 

Overall we wish to note the dramatic consequences of including a simple 
linear model of price feedback into a stochastic model. The outcomes of a hybrid 
of arithmetic and geometric Brownian motion are qualitatively more diverse 
than the possibilities predicted by cither in isolation, and a hybrid offers several 
features observed in the market. In particular variance explosion is predicted, 
as are power-law densities under certain conditions. Of course, much more 
realistic simulations might be done based on observed details of the order-book 
dynamics, better models of price-impact, and the endless non-Markov ideas of 
actual technical trading. However, while these might be readily simulated on a 
computer, the much-simplified approach presented here at least offers the hope 
of analytical tractability and insight. 
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